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Abstract. We give a sufficient condition on nonlinearities of an SDE on a compact connected 
, Riemannian manifold M which imphes that laws of all solutions converge weakly to the normal- 

■ ized Riemannian volume measure on M. This result is further applied to characterize invariant 

and ergodic measures for various SDEs on manifolds. 
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1. Introduction 



This work has its origins in an attempt (so far not concluded) by the authors to find a mathcmat- 
icahy sound numerical approximation for the stochastic geometric wave equations. Such equations 
have been recently introduced by two of us in [B], and our aim is to produce a counterpart of the 
1^ . numerical scheme introduced by three of us in [T] for the stochastic Landau-Lifshitz-Gilbert equa- 

te ' tions (used in the theory of ferromagnetism) . For the deterministic geometric wave equations such 

questions have recently been studied by two of us and Schatzle in [2] . During the initial stages we 
encountered problems related to large time behaviour of the solutions and we realised that we did 
not know an answer to such problems even for finite dimensional models. The reason being that 
the finite dimensional approximations to stochastic geometric wave equations, for instance those 
. generated by finite elements methods, are highly degenerate with respect to noise. Moreover, there 

CNJ ' are not many papers on the stochastic Langevin equations on manifolds, but see |16) and [3]. 

, On the other hand, the existence and uniqueness of invariant measures or convergence of so- 

lutions to an attracting measure was studied in case of uniformly elliptic diffusion operators e.g. 
, in [15] and in case of hypoelliptic parabolic and elliptic systems in 13 and [M]. The variety of 

' interesting problems is of course much wider than the cases described above as even quite simple 

SDE problems degenerate in such a way that none of the existing results covers them (see the 
examples at the end of this paper) . A particular approach to such degenerated elliptic problems 
is always required and we present one of them - still quite general - which implies that laws of all 
solutions converge weakly to the normalized Riemannian volume measure on M. For instance, in 
^ , Theorem 12.91 we formulate sufficient conditions for the existence and uniqueness of an attractive 

measure. We develop further this result to characterize invariant and ergodic measures for various 
degenerated SDEs on manifolds. To be more precise, in Theorems l3.4l and l3.6l we give a complete 
characterization of the set of invariant and invariant ergodic measures for certain classes of SDEs 
on the sphere and the special orthogonal group SO3. 

A further goal of this paper is to construct numerical schemes for solving certain classes of SDEs 
on manifolds, i.e., the finite dimensional stochastic Landau-Lifshitz-Gilbert (LLG) equation, and 
the geodesic equation on the sphere with stochastic forcing. The LLG equation has been widely 
studied in the physical literature, see for instance and references therein. They have also 
received some attention from the mathematical point of view, see [17]. General stochastic LLG's 
for non- uniform magnetisations have been studied by one of us and Goldys in [5], while Lelievre 
et al in [19] have described certain numerical schemes for SDEs with constraints. A convergent 
discretisation in space and time which is based on finite elements is proposed in [T; this scheme 
guarantees the sphere constraint to hold for approximate magnetisation processes, and thus inherits 
the Lyapunov structure of the problem. As a consequence, iterates may be shown to construct weak 
martingale solutions of the limiting equations. Main steps of this construction are detailed for the 
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finite dimensional LLG equation in Section [3.1l to then study long-time dynamics in Section [3.2l A 
corresponding numerical program for second order equations with stochastic forcing is detailed in 
Section |4.1[ which requires different tools; in particular, a discrete Lagrange multiplier is used in 
Algorithm B for iterates to inherit the sphere constraint in a discrete setting. Overall convergence 
of iterates is asserted in Theorem 14.31 which holds for this particular SDE on the sphere, but 
which may also be considered as a first step to numerically approximate the stochastic geometric 
wave equation. Again, computational examples are provided in Section to illustrate the results 
proved in this work, and motivate further analytical studies for computationally observed long-time 
behaviors which lack a sound analytical understanding at this stage. 

2. The problem and preliminary results 

Let M be a compact connected d-dimensional Riemannian manifold whose normalized Rie- 
mannian volume measure is denoted by A (i.e. X{M) = 1). Integration with respect to A will 
be denoted by dx. Let us assume that F and Gi, . . . , Gm are smooth vector fields on M. Let 
W^, . . . , W"^ be independent (^t)-Wiener processes on a complete stochastic basis {ft, J^, (^f ), P) 
and consider the following Stratonovich SDE 

m 

(2.1) du^ F{u)dt + ^Gk{u)o dW'' 

fc=i 

on M. 

The following (preliminary) results on SDE on manifolds in this section are generally known 
and can be found e.g. in [11], [15l or fT8|. Denote by P = (-Pt)j>Q the Feller semigroup associated 
with the SDE 1^^, i.e. 

Pth{x) ^ E{h{u'' (t))) , h£Bb{M), {t,x)eR+xM 

by Bb{M) we denote the Banach space of all real bounded and Borel measurable functions on M, 
and is a solution of (|2.ip such that 

w"(0) =x. 

Let us also denote by C(Af) the separable Banach space of all real continuous functions on Af, 
let us recall that the restriction of P to C{M) is a Co-semigroup on C(Af), denote by A its 
infinitesimal generator and recall the following characterization: 

C'{M) C D{A) 

(2.2) Af = Ff+-Y,Gk{Gkf), feC'{M). 

fe=i 

We will denote by P* the dual semigroup on the space Ai (M) of all R- valued Borel measures 
on Af, and by A* the infinitesimal generator of the dual semigroup. 

We convene that, throughout this paper, all measures on Af will be Borel, i.e. with the domain 
<^(Af). 

Let us also introduce the transition kernel 

(2.3) p(t, x,A)=F {{u''{t) eA}), ; t>0, xe Af , A e ^(Af ). 
Let I' be a probability measure on Af . We say that v is invariant provided 

(2.4) / p{t, x, A) dv = v{A), for aU t > 0, x G Af and A e ^(Af ). 
Jm 

An invariant probability measure is said to be ergodic if it is an extremal point of the convex set 
of all invariant probability measures on Af . 

A probability measure v on Af is called attractive if 

(2.5) \\mPtf{x)=( fdv, for ah a; € Af and / € C(Af ) 

Jm 



We will also say that v is invariant for the SDE l|2.H l that generates the semigroup P. 
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which, in other words, means that every solution u of (j2.ip converges in law to v irrespectively of 
the initial condition (even if u(0) is random). 

Let Lf) = {F,Gi,..., G™} and L„ = {[X, F] : X, F e UjCo ^jl for n(^n, and denote 

OO OO 

(2.6) ^ = span |J i„, = span (Lq \ {F}) U |J L„ 

n— n— 1 

where [X, Y] ~ XY — YX . Then we say that the vector fields F,Gi, . . . , Gm satisfy the hypothesis 
(H) if dim{Xp : X E ^} ~ d for every p G M (the Hormander condition), 
(F) if every / e C°°(M) is constant on M provided Lf = on Af for every L G ^p, 
(D) if div Gi = • • • = div Gm = div = on M, 

(C) if, for every I e N, there exists a finite dimensional subspace Ci of C'^{M) such that 
-^(G;) C Gi and spanlJ^gj^Gi is dense in G(Af) with respect to the supremal norm. 

Remark 2.1. Let us observe that the equation (|2.ip and the hypotheses (H), (F) and (C) do not 
depend on the Riemannian structure of the manifold M. In particular, we can, if necessary (and 
possible), change the metric on M in order the vector fields F,Gi, . . . , Gm satisfy the hypothesis 
(D). Of course, a change of the metric on Al changes accordingly the Riemannian volume measure 
on M. 

Remark 2.2. The hypothesis (F) is clearly satisfied if either the union of closures of the components 
of {p & M : dim{Xp : X G -Sfi?} = d} is a connected dense subset of M, or if any two points of M 
can be connected by a piecewise smooth curve consisting of integral curves of the vector fields in 
Jifp (see the examples in Section |3] and S]). 

Proposition 2.3. A probability measure v on M is invariant if and only if 

I Ahdu^Q, heC^{M). 
Jm 

Proof. See the identity (4.58) on p. 292 in [H]. □ 

Proposition 2.4 (Krylov-Bogolyubov) . There exists at least one invariant measure. 

Proof See Corollary 3.1.2 in [9]. □ 

Proposition 2.5. Let R G G°°{M) be a non-negative function. Then the measure dfj, ~ RdX is 
an invariant measure iff the following equality holds on M 

1 

(2.7) - div{[div(i?Gfc)]Gfe} = div(i?F). 

fe=i 

Proof. This result can be proved directly as for instance the last but one claim of Theorem 6.3.2 
from [TD] or it can be deduced from that result. □ 



Corollary 2.6. // (D) holds th 



en 



(2.8) A*f = -Ff+^J2^l(f^^ f&C'iM). 

k=l 

Moreover, the condition \2. 7p is equivalent to A* R = 0. 
Proof. Since for a vector field X on M and any / G G^(A/), we have 
X*f = -div{fX) = -fdWX-X{f), 
{X^Yf = f{dWXf + 2X{f)dWX + fX{dWX)+X^f), 
bv Proposition [2751 we get the result. □ 



Proposition 2.7. // (H) holds then every invariant measure has a density R G G°°{M) with 
respect to the Riemannian volume measure X on M . 

Proof. This follows from [T^]. □ 
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The following is a consequence of Theorem 6.3.2 from |10) (or of our Proposition 
Theorem 2.8. The normalized Riemannian volume measure X on M is invariant iff 

m 

(2.9) dw{F) + ^[Gk{dWGk)- {divGkf]^0 on M. 

fe=i 

In particular, if the condition (D) is satisfied, then A is invariant. Moreover, if both conditions 
(D) and (H) are satisfied, then A is the unique invariant probability measure. 

Proof. It is sufficient to prove the last claim. So let us assume that fi is an invariant probability 
measure. Then, in view of Proposition 12. 71 /i has a density R with respect to A and R e C°°{M). 
Moreover, by Corollarv l2.6[ A*R = 0. On the other hand, in view of the first formula in the proof 
of Corollary 12.61 from the assumption (D) we infer that 

(2.10) / R- A*Rdx = ^]-y^ j \GkR\'^dx. 

Hence GiR = • • • = GmR — and consequently FR = —A*R — 0. Thus 

ii? = for every L e if (i^, d, . . . , G„). 
The condition (H) yields that i? is a constant function. □ 

Theorem 2.9. Asuume the that the Hypothesis (F), (D) and (C) hold. Assume that a process u 
is a solution to SDE i2.1\) . Then there exists a probability measure 9 on M such that 



(2.11) lim E/(M(t)) = / fdO, feC{M). 

Jm 

If, in addition, (H) holds then 8 — X. 

Proof. Let us fix a natural number I S N. Let us choose a basis {/i, . . . , /„} of C; and let (aij)"^-^! 
be a a real {n x n)-matrix of the restriction of linear operator A to Ci with respect to this basis. 
In particular, 

n 

-4/j(p) = ^ay/,(p), p£M, i£{l,...,n}. 

i=i 

Let us denote by A the linear operator in M" whose matrix in the canonical basis {ei, • • • , e„} is 
equal to (ay )"j=i- 

Let a process u be a solution to SDE (|2.ip . Then for every i £ {1, . . . , n} 

E/,(M(t)) = Ef,{um+E f {Af,)iu{s))ds 

Efj{u{s)ds, t>0. 



n „t 

•-fMO))+Y.a,, / ] 



Hence 



(2.12) ^E/,(ti(i))e, = e^*[^E/j(ii(0))e,], t>0. 

i=i j=i 

Since by Lemma 12.101 below the operator valued function e"^* is convergent as t ^ oo, in view of 
(I2.12[) we infer that so are the functions E/i(u(t)). By the density part of assumption (C), we 
conclude that for every / £ C{M), Kf{u{t)) is convergent as t — oo. 

Now we will prove the last claim and so we assume that the condition (H) holds. Then by the 
Krylov-Bogolyubov Theorem in conjunction with Theorem 12. 8[ 

lim - / iys{B) ds = X{B), B £ 3§{M), 

t^ca t Jq 

where vt is the law of u{t). This concludes the proof. □ 
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Lemma 2.10. The operator valued function e*"^ from the proof of Theorem \2.9\ is convergent as 
t — >■ oo. 

Proof of Lemma \2.10\ Let us begin by introducing the following notation; 

n 

y{p) = {fj{v))U = E fM^i e P e M. 

If follows from equality (|2.12p that for any p £ M, supj>o |e'^*j/(p)| < oo. On the other hand, since 
the functions /i, . . . , /„ are linearly independent in C^(M), the vectors {y{p) : p G M} span M". 
Hence, we infer that 

(2.13) suplle-^*!! < oo. 

t>o 

In the above we considered A and e"^* as linear operators on R". We may naturally extend them 
to the complex space C". Since the euclidean norms of e"^* in M" and C" coincide, we infer that 

(2.14) sup||e^*||i(c.) < oo. 

It follows from (|2.14p that the spectrum a{A) is contained in C_, the closed real-negative halfplane 
of C. 

We will show that {iM. \ {0}) n cr(A) = 0. Arguing by contradiction let us consider A € 
(iK \ {0}) n o'{A). Then X = ia for some a e K \ {0} belongs to the spectrum of the adjoint 
operator A*. Hence and there is z = u + e C" \ {0} such that A*z = iaz, i.e. A*u = —av and 
A*v = au. If we denote g = J2j=i '^jfj ^^-^ h = '"jfj^ then A*g = —ah and A*h = ag. 

Hence, by identity (|2.10p . we get 

0=/" g[A*g + hA*h]dx = -l-^ [ [\Gjg\^ + \Gjh\^]dx. 

Therefore Gjg = Gjh = on M for every j e {1,...,™}. Consequently, Fg = —ah and 
Fh = ag and so 

[G,,F]g = G,Fg = -aGjh = 0, 

[Gj,F]h = GjFh = aGjg = 0. 

By an induction, we infer that Xg = Xh = for every X E ^p. Hence, by assumption (F), 
the functions g and h are constant on M. Thus h = —a^^Fg = and g — a~^Fh = what 
implies that z = 0. This contradiction concludes the proof of our claim. 

Let us now assume that S <j{A). Then there exists a natural number S N, matrices Ak, 
k £ {1, . . . ,N}, a positive number e > and a holomorphic £(M")-valued function V defined on 
an open disc in C of radius e such that 

N 

(A-^)-i -EA-'^ylfe + F(A), 0<|A|<e. 

k=l 

By the residuum theorem, 

k=l ^ ' 

where o(t) — > as t — > oo. We however already know that e*'^ is bounded and hence, necessarily, 
TV = 1. Thus we infer that e*"^ — > as i — )■ oo what concludes the proof. 

□ 
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3. Example I 

Assume that A, Bi, . . . , are antisymmetric {n x n)-matrices and Wi, . . . ,Wm are inde- 
pendent (^f)-Wiener processes. Consider the foUowing stochastic differential equation in the 
Stratonovich form, 

m 

(3.1) dz = Azdt + Y^ Bkz o dWk 

fc=i 

on the sphere §"^^ C M". Denoting by £(M") the space of real [n x n)-matrices, every solution of 
(|3.ip can be written as Z{t)zo, where Z is an £(R")-valued solution of the following Stratonovich 
stochastic differential equation 

(3.2) Z{t) ^ AZdt + J2T=iBkZo dW'', t>0 

satisfying the initial condition Z{0) = I. In fact, the solution Z takes values in the set SOn called 
the special orthogonal group and consisting of all unitary operators in R" with determinant 1. Let 
us recall that 

• the set SOn is a compact connected Lie group and a submanifold of £(R"), 

• for Z e SOn, 

TzSOn ^ {V e £(M") : V + ZV*Z = 0} 
= {y e £(R") : VZ* + = 0}, 

• there exists a bi-invariant Riemannian metric on SOn and the corresponding normalized 
Riemannian volume measure A on SOn is bi-invariant (with respect to multiplication by 
elements of SOn)- 

For an antisymmetric M e £(R"), we denote by M the following vector field on SOn, 

M : SOn 3 Z^ MZ e TSOn- 

Lemma 3.1. In the above framework, the conditions (D) and (C) are satisfied for the vector fields 
A and Bi, • • • , Bk- 

Proof. Let us fix an antisymmetric M £ £(R"). We have divM = 0. Indeed, if (Ei) are orthonor- 
mal vector fields on SOn, then, see [8j, 

(3) 

divM = ^(V£.M,£;,)R„x„. 

On the other hand, since M is antisymmetric, 

{\/EM:Ei}ur^.,^ =\iiRt'^{M{z + tE,z)- Mz,Eiz)jg,,..„ = {ME,z,E,z)v^^.r. =0. 

Combining the above two identities we infer that the condition (D) is satisfied. 

Let us prove the the hypothesis (C) is also satisfied. Towards this end we identify the space 
£(R") with R"^". Let us denote by C; the vector space spanned by the restriction to 0„ of 
polynomials on R"^" of order smaller or equal to /, i.e. Ci = spanj/o, : \a\ < I}, where 



n n 

nxn 



/„(X) = X" = []n<r' ^ = N]eR"^", aeW^ 

3 = 1 k=l 

Observe that 

n n 

(M/,)(Z)= "i,ia,fcz"-'^''+'^" = ^mkfc.-eo>'+e^''{Z), Z = [z,,]eW'', 

l,j,k=l l,j,k=l 

where e-''^ is the zero (n x r7,)-matrix except for the position (jfc) where e'j^ ~ 1. In particular, 
since the degree of fa-ei'^+e'^'' is < ^, we infer that Ci is invariant for the operator A defined in 
(|2.2p . Finally, by the Stone- Weierstrass Theorem, Ui^o dense in C{SOn)- □ 
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Corollary 3.2. The normalized Riemannian volume measure A on SOn is invariant with respect 
to the Feller semigroup generated by the SDE h3.2\) . 

Theorem 3.3. // the Lie algebra Ji'A{A, _Bi, . . . , Bm) contains all antisymmetric {n x n)-matrices, 
then A is an attractive and unique invariant probability measure for the Feller semigroup generated 
by the SDE ^Ol) . 

As a special case let us assume that A and B arc two antisymmetric (3 x 3)-matrices such that 
B 0. Consider the following Stratonovich SDE on SO3: 

(3.3) dZ{t) ^AZdt + BZo dW. 

Since both A and B are antisymmetric, we infer that 

(i) [A, J5] = iff j4 and B are linearly dependent, 

(ii) {A, B, [A, B]} spans the 3-dimensional space of antisymmetric matrices iff A and B are linearly 
independent 

Hence, two possibilities may arise. We begin with the first one. 

Theorem 3.4. // then the measure A on SO3 is an attractive and unique invariant 

probability measure for for the Feller semigroup generated by the SDE 113.3]) . 

If [A, — then A and B are commuting and so exp{tA + WtB} is the unique solution of 
([3:31) . Assume that B^O. Let us put p = {^WBWl^^^)^. For K e {S^, SO3}, define 

(3.4) S : X K 3 {p, Z) ^ s{p)Z e 

(3.5) s{x, y) = ^B' + ^B + I, (x, y) e S^, 

P P 

Note that for antisymmetric (3 x 3)-matrix B, B^ — —p^B. Thus the definition (|3.4II3.5[) of the 
map S is correct. Moreover, since s{p)s{q) — s{pq) for all p,q eS^, the map 

(3.6) 3p^ Sp = Sip, •) e K 

is a group homomorphism, i.e. S{p, S{q, Z)) — S{pq, Z) for p,q G and Z € K . For Z E K we 
will denote by S^ the orbit of the group (Sp) ^^gi through Z, i.e. the function 

(3.7) S^ -.S^ 3 p>-^ S{p,Z) e K 
Let also denote by H be the space of all such orbits, i.e. 

H^{S^[§^] -.ZeK}. 

The space H is equipped with the quotient topology for the surjective projection 

TT-.K 3 Z^ e H 

with respect to which H is compact. The quotient topology of H is metrizable by the classical 
Hausdorff metric 

p{X,Y) — max {sup d{x,Y), sup d{y,X)}. 

xex yeY 

Let us denote by dx the normalized Haar measure on S^. li X £ H and U,V £ K are two elements 
belonging to the orbit of X, then S^{dx) = S^ (dx) on ^{K). Therefore we can define a measure 
px S^ (dx) on K with the support on X. Moreover for every J e ^^{K), the function 

H 3 X^ px{J) e [0,1] 

is Borel measurable. 

If is a probability measure on K, then by = 7r(i/) we denote the image probability measure 
on H by the map tt. Finally, we define a measure P on X by 

v{J)^S{dx®v){J)^ ( px{J)dv.,{X), Je.^{K). 

J H 

In the following three theorems, the equation p.3|) will be considered on K e {S^,503}. The 
joint proof is given after Theorem 13.71 
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Theorem 3.5. Let us assume that [A, B] ~ and B =/= 0. Let Z be a solution of the equation 
p.3p on K Q {§^,5*03} and let us denote by v the law of Z{0). Then the law of Z(t) converges 
weakly as i —> 00 to the measure v on SS(K\ 

The following notation will be useful in formulation of the next result. If is a probability 
measure on i7, then vq is a probability measure on K defined via the following averaging formula. 

(3.8) MJ)= f PxiJ)de{X), Je.^iK). 

JH 

Theorem 3.6. Let us assume that [A^ B] = Q and B ^ 0. Then the mapping 9 t-^ ve is a bijection 
between the set of probability measures on H and the set of invariant probability measures for the 
equation (|3.3p on K Q {S'^jSO^}. 

Theorem 3.7. Let [A,B] — and B ^ 0. Then a probability measure v on K ^ {S^^SO^} is 
ergodic for the equation (|3.3|) on K iff there exists X ^ H such that v = [ix ■ 

Proof. Lets us observe that since B'^ ~ —pB, 

1 — cosf/os) „o sinffls) „ ^ 

(3.9) 6"^ = ^^5^ + — \!L±B + I, seR. 

P P 

Since A and B are linearly dependent, there exists n such that A = kB and. as already mentioned, 
every solution Z of the equation p. 31) has the form 

Z(<) = e*-^+'^'^Z(0) = e(«*+^')^Z(0) = S{-/{t), Z{0)) 

where 7(t) = (cos(Ki + Wt), sm{Kt + Wt))- Using the Fourier series, we may easily prove that ^{t) 
converges in law to the normalized uniform measure dx on ^(S^), hence Z{t) converges in law to 
i> = S{dx (E) where v is the law of Z{0) on 3i{K). Consequently, Theorem 13.51 is proved on one 
hand, and on the other hand, we proved that if is invariant then v ^ D. 

However, ® D) — v for every probability measure /i on ^(S^) (in particular, for the law of 
7(i)) hence u is invariant v = v. Moreover, S'(/i ® ve) — vg for every probability measure /i on 
^(S^) and every probability measure on £S{H). Hence {I'g : 9 probability measure on ^{H)} 
coincides with the set of invariant probability measures for the equation (j3.3l) . The injectivity of 
9 ^ vg follows from the fact that the image measure iT{vg) coincides with 9 on [^{H) and Theorem 
13.61 is thus proved. 

To prove Theorem 13.71 realize the following: If a probability measure 9 on M{H) is not a Dirac 
measure then there exists an open set Hq C H such that 9{Ho) G (0, 1). If we define Hi = H \ Hq 
and 9,{J) := 9{J f] H,)/9{H,) for i e {0, 1} and J G ^{H) then vg = 9{Ho)yg„ + 9{Hi)ue^ and 
thus vg is not extremal. On the other hand, if = ^Sx — ^o^Oo + '^i^Oi — ^Xoeo+MBi holds 
for some X G H, > 0, Xq + Xi — 1 and for some probability measures 6*0, 9i on ^{H) then 
Sx = '^i^Sx) ~ ■"'(^Aoeo+Aiei) = '^o^o + ^^i^'ii hence 9o — 9i and iix is thus extremal. □ 

Remark 3.8. If [A, B] ~ 0, B Q and is a probability measure on K then the following three 
conditions are equivalent. 

(i) v is an invariant measure for the equation (j3.3p . 

(ii) Sp{h') = v for every p G 

(iii) 1/ = S{dx (g) v). 

Remark 3.9. In a finer look, Theorem 13.51 tells us also what happens, as t 00, on fibres of tt. 
If -B] = then the diffusion p.3p uniformizes, as t — ?► 00, the initial distribution v along the 
orbits = {S{p, Z) : p € S^} indexed hy Z G K. Indeed, let j/* be the law of the solution (with 
the initial law at time t > and consider the j^*-unique system of conditional probabilities 
(j/* [.[tt = X])xeH = {i^x)xeH on ^{K), i^^iX) — 1, to which the the measure desintegrates 
with respect to v^. Then it can be verified by the definition of a desintegrated measure that 
— S{lt (8) I'x) where k is the law of ^{t), hence converges weakly to fix by Theorem 13.51 
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3.1. Numerical approximation. The stochastic Landau-Lifshitz-Gilbert equation describes (uni- 
form) atomistic ferromagnetic spin dynamics at finite temperatures; it is of the form p.3p . with 
Az — —z X h and Bz = —z x [h + h±), such that 

(3.10) dz = -z X hdt- z X {h + hj_)dW , t>Q, 

with z{0) — zq. Here, ZQ,h e M.^ are some unit vectors, and h± e M.^ be perpendicular to h. The 
Hamiltonian £(ip) = —{h,(p)^3 is conserved by the flow for absent stochastic forcing, or h± = 0. 
We propose a non-dissipative, symmetric, and convergent discretization of (|3.10|) . Let Ik be an 
equi-distant mesh of size A: > covering [0, T]. We denote :— ^{(p"~^^ + (p"). 

Algorithm A. Let Z° zq. For neN, find e ^^^^ ^^^^t 

(3.11) - Z" = -fcZ"+i/2 xh- Z"+i/2 X (/i + /i^)AVK„+i . 
where AW^+i := t^(Wi) " W^(^n) ^ -^(0, fc)- 

Let fc < 1. Similar to [1], the R^-valued random variables {Z"'^^}n>o exist, are unique, satisfy 
|^n+i| _ 2 for all n > 0, and converge to the solution of p.lOp for k 0. Moreover, pathwise 
conservation of energy f (Z"+^) = £{zq) holds for h± = 0. 

Remark 3.10. Let := EZ''+\ Then for 7^ := jHI^, 



(3.12) ir"+i - jT" + /cjr"+i/2 X /i + + /i_Lp(^"+^/^ - (ir"+l/^ 



Hence, the stochastic forcing exerts a damping in direction perpendicular to h+h±. Estimate (j3.12p 
follows as in the (more detailed) Remark 14.41 below for a second order equation with stochastic 
forcing. The main ideas are to repeatedly use p. lip , and approximation arguments that base on 
= 1 As a consequence, the limiting equation reads 

3^t + ^><h+^\h + h^\^{S'-{S',h)h) =0. 

3.2. Numerical experiments. We perform computational studies of the stochastic Landau- 
Lifshitz-Gilbert equation p.lOp using Algorithm A. The nonlinear system in Algorithm A is solved 
up to machine accuracy by a fixed-point algorithm, cf. 1 . By N we denote the number of discrete 
sample paths of the Wiener process used in the computations, and denotes the numerical 
solution at time computed for the l-th sample path. 

The unit sphere is divided into segments ujij C associated with points 

x.y = (sin(i7r/16) cos(j'7r/16), sin(i7r/16) sin(j7r/16), cos(i7r/16)) , 

i — 0, . . . , 16, J = 0, . . . , 31 such that ujij = {x G E>^\ = argminx,^ |x — x/m|}. For the above 
partition of the sphere, at a fixed time level i", we construct a piecewise constant empirical 
probability density function /"(x) : — > i? as 



r(x)u., -r(x.,) 



for i = 0, . . . , 16, j = 0, . . . , 31, where #i7 denotes the cardinality of the set fi. 

The presented results below were computed with k = 0.01 for N ~ 20000 sample paths. 

For the first numerical experiment we consider zq = (0, 1/a/(2), l/-\/(2)), h — (0,0,1), h± — 
(0,0,0), which corresponds to [A, i?] = 0. The resulting probability density for any zq is uniform 
on a circle with the center {zo,h)h, see Figure ??. Further, E(z(<)) — ^ {zo,h)h for t — >■ oo and 
{z{t), h) = (zq, h) for any zq (i.e. the pathwise angle between z and h is constant in time) for any 
t. The computations agree with the statements of Theorems 13 . 511X71 

In the second experiment we set zq — (0, l/\/2, l/\/2), h = (0,0,1), h± = (0,1,0). In this 
case we have [A, B] ^ and the initial probability density which is a Dirac measure concentrated 
at Zq is expected to approach the uniform probability density function of the unit sphere — 
(47r)~^ « 0.0796. We also expect that E(z(t)) — >■ for t ^ oo. Because of the finite approximation 
of the problem, the uniform probability density is only attained approximately. Once the 
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computed probability density becomes diffuse, it fluctuates randomly around the uniform state 
/■^j see Figure ??. Analogically, the trajectory E{u{t)) approaches the center of the sphere and 
for long times fluctuates randomly around the center, see Figure ??. The random fluctuations 
in the probability density function /" can be significantly reduced by taking the average over a 
sufficiently large number of time levels. Here, we compute the time averaged probability density 
function / over the last 100 time levels, i.e., we take /(x) = -j^ ST/fe-ioo /"(^)' Figure ??. 
The results in the second numerical experiment agree with the assertion of Theorem 13.41 

4. Example II 

The geodesic equation on the sphere has the form dii = —\u\'^udt. Consider this second order 
equation with a stochastic perturbation 

(4.1) dii = -\ufudt + {u X ii) o dW 

which is a Stratonovich SDE dz = F{z) dt + G{z) o dW on the tangent bundle TS^ C driven 
by a standard one-dimensional Wiener process W where u x li e M'^ is the outer product in M'^ 
and 

are tangent vector fields on TS^. Observe that F, G and [G, F] are also orthogonal tangent vector 
fields on the 3-dimensional submanifold 

(4.2) Af, = {(p,e)er§':|ei=r}, r>0 

hence the hypothesis (H) is verified. If we put Ei — F/{r^ + r^)^, E2 — G/r, E3 — [G,F]/r then 
{El, E2, E3} is an othonormal frame on Mr and 

3 

divr = ^(dfi^.r, Ej)^e =0, Ye {f, g, [f., g]} 

where dxY{p) = \imt^ot^^[Y{p + tX) — Y{p)] and so the hypothesis (D) is holds. Next, if 
[G,F]f — Gf = on Mr for some / G C°°{Mr) then / is constant along any integral curve 
7 = [G, F]j or 6^ GS. Since 

Rng7 = {(p,0eM, :p = 7i(0)} and Rng <5 = {(p, G : ^ - (52(0)}, 

any two points in Mr can be connected by a piecewise-smooth curve consisting of at most three 
integral curves of [G, F] and G hence the hypothesis (F) is satisfied. 

Finally, let Ci denote the space spanned by polynomials on of order smaller or equal to I. 
Then F[Ci] U G[Ci] C Ci hence Ci is invariant for the operator A defined in (|2.2I) and Uj^o 
dense in C{Mr) by the Stone- Weierstrass theorem and the hypothesis (C) holds. 

Let Hr be the probability measure on TEP supported on Mr for r > such that the restriction 
of fj,r on ^M(Mr) coincides with the normalized Riemannian volume measure on Mr- If is a 
probability measure on T§^, we define the measures 

MU) ^iy{ip,0^ : 1^1 e C/}, U (E ^{0, 00) 

D{A) = i^{AnMo)+ I pr{A)du^, Ae^(T§2). 
J (0,00) 

Theorem 4.1. Let {u,u) be a solution of j| ) on T§^ with an initial distribution v on SS{T^\ 
Then the laws of (m, u) converge weakly to 9 as t ^ 00. Moreover, v is invariant for ( |^._?| ) iff 
V — V and v is ergodic for j[ ) iff v = 5(pfl) for some p G §^ or v — fir for some r > 0. 

Proof. The spaces {(p, 0)}pgs2 and Mr for r > are invariant for (j4.ip . Hence, if / e G(T§^) 
then f dp {t,x) converges to f{p,0) if a; = (p, 0) or it converges to Jrp^2 f dfir if a; = {p,£,) 
and 1^1 = r > by Theorem 12.91 since fXr is the unique invariant and also attractive probability 
measure for (|4.ip on Mr. The ergodicity follows from Remark 14.21 as ergodic probability measures 
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are the extremal points of the set of all invariant probability measures (see e.g. Proposition 3.2.7 
in 0). □ 

Remark 4.2. Invariant measures for ()4.ip can be uniquely described as measures 



ly^^hiA) = a n Mo) + / firiA) db, A e ^{T§^) 

J{0,oo) 

where a and b are finite measures on .^(Afg) and on ^(0, oo) so that a(Mo) + 6(0, oo) ~ 1. 

4.1. Numerical approximation. We propose a non-dissipative, symmetric discretization of 
()4.1|) to construct strong solutions and numerically study long-time asymptotics. Let {{U''\ V")}n 
be approximate iterates of {{u{tn),u(tn))}n on an equi-distant mesh 1^ of size fc > 0, covering 
[0,T]. We denote dtip"+^ := - v")- 

Algorithm B. Let iU°,V°) (uq, u(0)) , and ~ f7"-fcV°. For n>0, find {U"+\V"+'^ , X" 
R3+3+i^ such that for AW^+i W{t„+i) - W{t^) - ^(0,fc) holds 



\n+l 1 

+ 1 



^n+l_-j^n ^ ^ (J7"+l + C/"-l) + -(C/"+l + [/"-l) X (y"+l+F")AVF„ 



2 ' '4 
(4.3) dtU"+^ = 

for i([/"+i + [/"-!) = 0, 

A"+i = <! for i(C/"+i + [/"-!) ^0 and ,i> 1, 

for i(C/i + [7-1)^0 and n = 0. 

The choice of the Lagrange multiplier A"+^ ensures that |?7"+^| = 1 for n > Q; the case n = has 
to compensate for the fact that the defined U"^ is not necessarily of unit length. 

To see the formula (gSlla for n > 1, we muhiply (|431) with i([/"+i + C/""^) and use the discrete 
product formula 

{dtcp\ r) = ~{v''-\dtr) + dt(^", r) 

to find 

where (•, •) denotes the scalar product in R^, and | • | = (•, ■f'^. Since (\/"+\ C/"+i + f/") = 0, we 
further obtain 

= -i(^(F", V""+^ + F""^) - (F"+\ F") + (V^",F""^)^ = -(F", F""^) . 

Hence -(F", V"^) = A"+i|i(t/"+i + C/"-i)|2, which yields the formula for A"+i in 
For n = 0, we conclude similarly, using (f7°, F*^) =0, and the definition of C/^^. 

Theorem 4.3. Let T > 0, and fc < kf){U^, V*^) be sufficiently small. For every n > 0, there exist 
unique MP^^-valued random variables ([/"+^, 1^"^^) of Algorithm A such that \U"^^\ — 1 for all 
n > 0, and 

E{V''+^) = E{V°) where E{ip) = . 

Moreover, {{U"^^ ,V"'^^)}n>o construct strong solutions of for fc — > 0, in the way specified 

in the proof below. 



Solvability of ()4.3p will be shown by an inductional argument that is based on Brouwer's fixed 
point theorem: an auxiliary problem is introduced which excludes the case where i(?7""''^ + 
f/""^) = when computing A"+^; for sufficiently small fc > 0, constructed solutions of the 
auxiliary problem are in fact solutions of (|4.3p . 
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Proof. 1. Auxiliary problem. Fix n > 1. For every < e < |;, define the continuous function 
: _^ R3 wliere 

(4.4) F^{W) ■.= j{W- [/") + k ^^"'i!^" -Wx (F" - ^;7")AW^„+i . 

We compute respectively, 

|(w^-;7",w^) > l{\w\^\u^\)\w\, 

lyn yn-l\ 

max{|M^|^, ej 



Since the stochastic term in (j4.4l) vanishes after multiplication with W , there exists some function 
i?„ > such that 

Then, Brouwer's fixed point theorem implies the existence of W* = ^(C/"+^ + C/"^^), such that 
J","(i([/"+i + = for aU u} efl. 

The argument easily adopts to the case n = 0. 

Solvability and energy identity. We show that i([/"+i + J7"~i) solves {^{U"+^ + U"-^)) = 
0. By induction, it suffices to verify that 

|i(c/"+i + c/"-i)| = |-y"+i + i(f7" + [/"-i)| >!-([/" + 

2 2 2 2 2 

(4.5) > l-\-'l\V-+'\>l, 

for fc < ko{U'^,V'^) < 1 sufficiently small. 

Let n>l. For aU < ^ < n, there holds |C/^| = 1, and 

(4.6) EiV^) = S(yO) . 
Then W* = i(C/"+i + C/"-i) from Step 1. solves 

(4.7) fcdty"+i = ([/"+' + C^""') + + U''-^) X (!/"+! + V'')AWn+i , 
where 

A"+i = - 



(T/»,F"-i) 
max{e, |i(f/"+i + C/"-i)|2} 



Testing gH]) with i(J7"+i - C/"-i) = |(V^"+i + F"), and using binomial formula i(C/"+^ + 
[/"-!, C/"+i-J7"-i) = i(|C/"+i|2-l), as weU as < fc2|T/"+i|2 + |J7"|2, and the inductive 

assumption = 1, 

rf,|l/"+l|2 < ^^fc2|^n+l|2 

|T/n| |T/n-l| i.2|T/n+l|2 1.2 

4max{e, (1 — I — fc|y"+^|)2} 4e 

for e < i. By a (repeated) use of the discrete version of Gronwall's inequality, there exists a 
constant C > independent on time T > 0, such that 

(4.9) < C|V^°|2. 

As a consequence, g3]) is valid, and hence J","^(i([/"+i + C/"-i)) = for e = 0; therefore, C/"+i 
solves (14. 3p . satisfies the sphere constraint, and conserves the Hamiltonian, i.e., (|4.6p holds for 
0<i<n + l. 
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For n = 0, we argue correspondingly, starting in (j4.8l) with 



max{e,|i(C/i + C/-i)|} 



from which we again infer (|4.9p , and ()4.5p . 

Uniqueness of {(t/", follows by an energy argument, using ()4.5p . (|4.6p . k < ko, and the 

discrete version of Gronwall's inequality. 

5. Convergence. We rewrite (j4.3p in the form 

dv = —\v\'^udt + uxvodW, 

(4.10) du = vdt, 

For K-^-valued iterates {</3"}„>o on the mesh Ik that covers [0,T] we define for every t S [i„,i„+i) 
functions 

^feW V^", and ^+(0 

We now show that corresponding functions of Algorithm A satisfy P-almost surely 

(4.11) (^±,r±) ^ (w,w) inC([0,T];R) (fc' -> 0) , 
where {u,v) is strong solution of ()4.3p . 



It is because of the discrete sphere constraint and the (discrete) energy identity that sequences 

{(^fe^n^)},cC([0,T];R) 
are uniformly bounded. Moreover, there holds for all t > 

/■*+ \ + 

y{t) = r(o)+ / — [^+ + '^--fcr-]ds 

Jo 2 

1 r'^ 

(4.12) +- / ('^+ + -ky-)x {■r+ + y-)dw{s), 



4 .„ 

•^(t) = ^(0)+ / r(s)ds. 
Jo 

Then, by (14.12^ 9. (j4.5p . and Holder continuity property of VF, sequences {('^fc,'^)}^, are equi- 
continuous. Hence, by Arzela-Ascoli theorem, there exist sub-sequences {{^k' ,%:')} ^ ^^d con- 
tinuous processes ('^, f) on [0,T] such that 

(4.13) ||^fc'-u||c([o,T];R3) + ||rfc'-w||c([o,T];M3)^0 (fc'^O) P-a.s. 

We identify limits in (|4.12l) . The only crucial term is the stochastic (Ito) integral term which may 
be stated in the form 

(4.14) -/ {'^+ + '^~ - ky-) X {y- + -f')dw{s) . 

2 Jo 2 
We easily find for every t e [0, T], 

1 f*^ /■* 

-/ (^+ + '2'--A:r-)xr-dW^(s)^ / uX'ydW^(s) (fc ^ 0) P-a.s. 

2 Jo Jo 

The remaining term in ()4.14p involves ^Y, which will be substituted by identity (|4.3P i . 

-J {'^+ + '^- -kr-)x [y- + k—{9/+ + '^- -ky-) 

(4.15) -f^^" - /c^-) X {y- + ^y)AW„+i^ dW{s) . 
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If compared to (I4.14p . the critical factor -l"^ is now scaled by an additional AWn+i', using again 
(14.3^ 1 ■ Ito's isometry, and the estimate EjAWn+ip" < Cfc^" ^ then lead to 

i ^ (<^+ + '^-)x ((<^+ + - ky-) X (r- + ^r)) {w+ - w-)dw{s) 



(4.16) 



u X {u X v)ds {k 0) 



for all t e [0,r]. As a consequence, there holds 



v{t)=v{0)+ / \v\^uds+ / uxwdW^(s) + 



1 



u X (u X v) ds , 



where the last term is the Stratonovich correction. The proof is complete. 

Remark 4.4. 1. Let \V°\ be constant, and := (EF"+\ EC/"+i). Then 



□ 



(4.17) 



< Ck\ 



i.e., the stochastic forcing term exerts a damping in direction To show this result, we start 

with 

(4.18) =: / + //. 

We use Theorem 14.31 and an approximation argument to conclude that 

/ = --Ef(y" F"-Mfi - fi 1 + [/""^)] 

= -^E[(y",i/"-i)('^"+i + -^"-1)] +c'(fc3) 



= -/sE|y"|^'2^ 



0(fc2), 



thanks to the power property of expectations, and |||(f/"'^"'^ + ?7"~^)P — l| < Ck^. 

We use the identity [/"+^ = kV^~^^ + If^ for the leading term in //, and properties of the vector 
product to conclude that 

// = ^E[(l/"+i - X V'AWn+i] + Je[([/" + U''-^) X (y"+i - y")AW„+i] =: Ih + Ih ■ 

We easily verify \IIa\ < Cfc^, thanks to (|4.3P i . and properties of iterates given in Theorem 
For //{,, we use (I4.3l) i as well, and the relevant term is then 



1e 

16 



(f/" + U''-^) X (([/"+! + u''-^) X (y"+i + y")) I A„+i 1^ 



-E 
2 



t/" X ([/" X y")|A„+i|2 +0(^2) 
= _ V" + ©(/c^) = + 0(A:2) ^ 

thanks to the power property of expectations, earlier boundedness results of iterates {(C/", V^")}ti, 
the fact that |J7"| = 1 for all n > 0, the cross product formula a x [b x c) = b{a, c) — c(a, 6), and 
another approximation argument. This observation then settles (j4.17p . 
2. Strong solutions of (|4.1I) satisfy 

|7/(i)| = i, E{v{t)) ^ E{vo) vte[o,r]. 



and are unique, due to Lipschitz continuity of coefficients in \A.1\ : hence, the whole sequence 
{('%, '^•)}fc converges to (w, u), for fc — > 0. 

3. Increments of a Wiener process in Algorithm A may be approximated by a sequence of 
general, not necessarily Gaussian random variables, which properly approximate higher moments 
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of AW,i_|_i; martingale solutions of (I4.ip may then be obtained by a more involved argumentation 
from Algorithm A as well, using theorems of Prohorov and Skorokhod; cf. [3] • 

4.2. Numerical experiments. In this section we present some numerical obtained by the Algo- 
rithm B that has been applied to a slightly more general problem than ()4.1|) 

dii = ~\u\'^udt + ^/15{u X u) o dW , 

where Z? is a fixed constant that controls the intensity of the noise term. The Lagrange multiplier 
was computed as 

The above formula is equivalent to the corresponding expression in (|4.3p . However, the present 
formulation (j4.19l) is slightly more convenient for numerical computations, since it ensures that 
the round off errors in the constraint |C/"| = 1 do not accumulate over time. The solution of the 
nonlinear scheme (j4.3l) is obtained up to machine accuracy by a simple fixed-point algorithm, cf. 

The probability density function /" was constructed as in Example I with N — 20000 sample 
paths. For all computations in this section we take the time step size k — 0.001 and the initial 
conditions u{0) = (0, 1, 0), u{0) — (1, 0, 0). The initial probability density function associated with 
the above initial conditions is a Dirac delta function concentrated around u(0). 

In Figure ?? we display the computed probability density /" for = 1, T = 60 at different time 
levels. Initially the probability density function is advected in the direction of the initial velocity 
and is simultaneously being diffused. For early times, the diffusion seems to act predominantly 
in the direction perpendicular to the initial velocity. In Figure ?? we display the time averaged 
probability density function /, the trajectory E(u(t)) and a zoom at ¥,(u(t)) near the center of the 
sphere. 

The evolution of the probability density for D = 10, T = 60 is shown in Figure ??. Similarly 
as in the previous experiment the probability density function diffuses and becomes uniform for 
large time. Some advection in the direction of the initial velocity can still be observed, however, 
the overall process has a predominantly diffuse character. We observe that the overall evolution 
damped due to the effects of the random forcing term, see Remark 14.41 (1) and Figure [T] In 
Figure ?? we display the time averaged probability density function, the trajectory K{u{t)) and a 
zoom at E{u{t)) near the center. 

Figure ?? contains the computed trajectories of E(u(t)) for D = 0.1 and D = 100. The 
respective probability densities asymptotically converge towards the uniform distribution for large 
time. 

In Figure [T] we show the graphs of the time evolution of the approximate error f^^x ■ ^" ~^ 
maxxgs2 |/"(x) — /■^l for D = 0.01, 0.1, 1, 10, 100. The quantity iS'/^ax serves as an approximation of 
the distance from the uniform probability distribution in the norm. Note that the oscillations 
in the error graphs are due to the approximation of the probability density. The numerical 
experiments provide evidence that the probability densities for all D converges towards the uniform 
probability density for t — > oo. The probability density evolutions for decreasing values of 
D have an increasingly "advective" character and the evolutions for increasing values have an 
increasingly "diffusive" character. It is also interesting to note, that the convergence towards the 
uniform distribution becomes slower for both increasing and decreasing values of D. 



Figure 1. Convergence of the probability distribution of u towards a uniform 
distribution for different values of the coefficient D. 

In the last experiment we study the long time behavior of the pair {u,u) for _D = 1, = 
20000. Towards this end, we introduce a partition of the manifold Mi defined in (|4.2|) . First, 
we consider a partition of the unit sphere into segments wf, i = 1,...,6 associated with the 
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points — (±1, 0, 0), (0, ±1, 0), (0, 0, ±1) in such a way that x G §^ belongs to ujf if and only if 
\x — xf\ — mini<j<6 \x — Xj\. Next, we denote by T; the tangent planes to points xf . Fixing an 
i G {1, . . . , 6}, the orthogonal projections of vectors {xf , . . . , a;f } onto the tangent plane Ti delimit 
4 sectors on T^. We subsequently halve each sector obtaining thus 8 equi-angular sectors 7,^^, . . . , 7^^ 
in Ti. Now we introduce the following partition of Mi into 6x8 segments (see Figure ??): a point 
(p,^) S TEP belongs to Mf if p G ujf and the orthogonal projection of ^ onto the tangent plane 
Ti belongs to the sector 7^ . It can be verified by symmetries of this partition that the normalized 
surface volume of each M/ is equal to 1/48. For n ~ 60000 (i.e., at time <" = 60) we have for 
i = 1,...,6, j = 1,...,8 #{l\Up G wf} G (3380,3260) w 7V/6 = 3333 and #{;| ([/;", V;") G M/} G 
(386,455) « N/6/8 ~ 417, see Figure ?? left and Figure ?? right, respectively. The numerical 
experiments indicate that the point-wise probability measure for (u, ii) converges to the invariant 
measure V. The (rescaled) approximate L^o error fmax for (u, u) has similar evolution as the 
approximate Loo error for u. Moreover, it seems that the convergence of the error in time is 
exponential, see Figure ??. 

Appendix A. Lie algebra 

Let U be an open set on a C°°-manifold. 

• The set ^ of all smooth tangent vector fields on [/ is a vector space with the Jacobi 
bracket. Any vector subspace of ^ closed under the Jacobi bracket is called a Lie 
algebra. 

• If A" is a set of smooth tangent vector fields on U, then we denote by ^{X) the smallest 
Lie algebra containing X. 

• U Ac ^ and p eU, then we define A{p) = {Ap : A G A}. 

Proposition A.l. Define Lq ~ spanjA"} and L„ ~ span{L„_i U {[A, _B] : A,B E i„_i}}. Then 
Proposition A.2. Assume that X C^. Let Xi, . . .,X,n,Y G ^ and let U G C°°([/). Then 

771 

^(Ai,...,A„,y)(p) = if(Ai,...,x„,y + ^/,x,)(p), peu. 

Proof. Let us write A^^{Xi,..., Y}, A^ ^ {X^, . . . , X^,Y + /, A^}, 

= 1^ ''"'^^ • ^ C°°iU), Lk G .^{A'), G n| . 

Apparently, is a Lie algebra for i G {1, 2}, A' C whenever {ij} = {1, 2} hence ^{A') C 
whenever = {1,2}. But then 

^{A')ip)c<^^{p)c^{Anip)- 

□ 
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